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SUMMARY 


The technique of modulation, or variable coefficients, is discussed 
and the analytical formulation is reviewed. Representative numerical 
results of the use of modulation are shown for the lifting and nonlifting 
cases. These results include the effects of modulation on peak accel- 
eration, entry corridor, and heat absorption. Results are given for 
entry at satellite speed and escape speed. The indications are that 
coefficient modulation on a vehicle with good lifting capability offers 
the possibility of sizable loading reductions or, alternatively, wider 
corridors; thus, steep entries become practical from the loading stand- 
point. The amount of steepness depends on the acceptable heating pen- 
alty. The price of sizable fractions of the possible gains does not 
appear to be excessive. 


INTRODUCTION 


For low entry angles, provision of a small L/D is helpful in 
the reduction of peak loadings. At the higher entry angles, defined 
arbitrarily as those greater than 6°, even high L/D will not keep 
the peak loadings below some nominal limit such as lOg. For super - 
circular entry a g limitation defines a corridor width. To stay within 
a g limit at steep entry angles, the technique of modulation, or vari- 
able coefficients, must be introduced. In this paper numerical results 
will be shown only for entry at satellite and escape speed; however, 
the indicated trends are independent of the velocity. 

In entry from escape speed, as in lunar return, the ability to 
withstand steep entry provides a margin for error in the guidance prob- 
lem; however, at satellite speeds large retro-rockets are needed to 
achieve steep entry angles . From a practical standpoint there are 
certain military purposes, emergency conditions, and abort conditions 
for which it is possible that steep entry will be required at satellite 
speeds . 



SYMBOLS 


drag coefficient 

lift coefficient 

re suit ant -force coefficient 

drag 

acceleration due to gravity 
altitude 


lift 

total heat absorbed 
ballistic modulation ratio 

lifting modulation ratio 

radius from center of earth to atmospheric entry point 
velocity 

horizontal distance 
angle of attack 

atmospheric scale height, - — 

d log e p 

flight-path angle with local horizontal, positive up 
minimum acceleration index 
parameter (ref. l) 


air density 
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Subscripts 

j 

MAX 

maximum 

MIN 

minimum 

Mod 

modulated 

o 

initial conditions 

Opt 

optimum 

pe 

perigee 

Unmod 

unmodulated 

1 

start of modulation 

2 

end of modulation 


A dot over a symbol indicates a derivative with respect to time. 
Bar over symbol indicates values for ( l / d )mAX‘ 

DISCUSSION 


The modulation principle is illustrated in figure 1 which shows a 
loading history with and without modulation for a lifting vehicle. The 
resultant aerodynamic force in units of the vehicle weight is plotted 
against the corresponding times. In the upper right corner, the corre- 
sponding trajectories are sketched. For the unmodulated case, the 
loading rises to a sharp peak just before the bottom of the pull-up at 
which point the curve has been ended (corresponds to the upper trajectory 
in the sketch) . For the modulated case, at some point in the entry the 
loading is not allowed to increase further and is held constant by con- 
tinuous reduction of the resultant -force coefficient (corresponds to the 
lower trajectory in the sketch) . Deeper penetrations of the atmosphere 
always occur in modulated pull-ups. If the proper point for start of 
the modulation has been selected, the vehicle will have nearly zero 
lift coefficient as it levels out. If the modulation is started too 
soon, the vehicle will not have leveled out sufficiently as the lift 
coefficient approaches zero. Higher loadings will be experienced sub- 
sequently than were maintained during the modulation period. For late 
starts of the modulation and consequent higher loading levels, the lift 
coefficient will not have approached zero at the bottom of the pull-up 
and the full capability of the vehicle will not be realized. For the 
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case shown, the shaded area next to the unmodulated curve may be loosely 
regarded as an impulse which is shifted to later times and lower loading 
levels for the modulated case. For precision, components of the impulse 
must be considered. The basic idea, however, is to replace a large force 
acting for a short time by a smaller force acting for a longer time to 
make essentially the same change in velocity. ri he modulation scheme 
indicated in figure 1 is that used in reference:-, 1 and 2. 

The following simple relation governs the air loading: 


Resultant force 
Weight 


« C R pV 2 


As indicated in figure 1 the rise in pV^ during entry is combated during 
the modulation period by the decrease in C R . The following two cases can 

be distinguished: 


Ballistic: 


C R - C D 


Lifting: 

c R = fz 2 + c L 2 

In the ballistic case the drag coefficient is reduced during modulation. 
In the more complicated lifting case both lift and drag coefficients are 
changed simultaneously according to whatever functional relation exists 
between them. 


The nature of the difference between the tvo cases can be visualized 
more concretely in terms of force polars. Force polars for the two cases 
are sketched in figure 2 along with the physical means necessary to 
attain the polars. In the ballistic case an axi symmetric vehicle is 
visualized as decreasing its frontal area from £ large value to a small 
value. The corresponding polar is part of the C D axis. In the lifting 

case a winged vehicle is visualized as changing attitude from an angle of 
attack of 90 ° to an angle of attack of 0° while the force coefficients 
trace out a loop in the Cj),Cl plane. This is the extreme range of 
lift modulation; the optimum range lies, as will be shown, between maxi- 
mum lift coefficient and minimum drag coefficiert. 

In the ballistic case a geometry modulation is needed; in the lifting 
case an attitude modulation is needed. Other cases can be imagined in 
which both geometry and attitude modulation occur simultaneously and thus 
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trace out paths intermediate to the two shown. It can he shown (refs. 3 
and 4) that the g-alleviation capability for the lifting case depends 
strongly on a single parameter, the maximum lift coefficient in units 
of the minimum drag coefficient. For the ballistic case, the allevia- 
tion depends on the ratio of maximum drag coefficient to minimum drag 
coefficient. (See ref. 1.) For unit minimum drag coefficient the 
parameters appropriate to the two cases are marked on the polars as r B , 

the ballistic modulation ratio, and r^, the lifting modulation ratio. 

In the analysis of reference 2, the modulation was effectively restricted 
to the portion of the polar on the low -drag side of maximum lift-drag 
ratio (between the two lowest dots) . This restriction leads to an under- 
estimation of the possible g alleviation and an overestimation of the 
heat load associated with a given g level. 


ANALYTICAL FORMULATION 


Isolation of the two modulation parameters is made possible by the 
separation of the modulation problem into two limiting cases with appro- 
priate analytic assumptions. (See fig. 3 .) Formulas derived with these 
assumptions are presented in appendix A. In the ballistic case varia- 
tions in the flight -path angle are presumed to be unimportant whereas 
the velocity changes are considered to be significant. In the lifting 
case the complementary assumptions are made; that is, changes in velocity 
are ignored whereas changes in flight-path angle are considered to be 
important. In physical terms, the assumptions correspond to separation 
of the two possible operations on the velocity vector, changes of magni- 
tude and changes of direction. On course, these assumptions do not apply 
to the portions of the trajectories beyond peak g which are indicated in 
figure 3 as dashed lines. 

The assumption of constant velocity in the lifting case is conserva- 
tive in the sense that it leads to an underestimation of the g allevia- 
tion. Neglect of the actual velocity decrease means that in the real 
lifting case higher lift coefficients can be used at the same altitudes 
for the same loading limit. A lower loading limit than is indicated by 
the analysis is thus possible in the real case. 

The ballistic case has been analyzed in reference 1 and the princi- 
pal result is indicated in figure 4. Plotted vertically is the peak 
resultant force in units of the peak resultant force for no modulation. 
The sole parameter on which the g alleviation depends is the modulation 
ratio r£, the ratio of maximum drag to minimum drag. A region of 

rapid initial gains is followed by one in which diminishing returns have 
set in. The alleviation ratio always decreases but always at a slower 
rate. The ballistic case is completely solved in terms of the ballistic 
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modulation ratio. In the lifting case, a different modulation ratio, 
maximum lift coefficient to minimum drag coefficient, plays the leading 
role . 


The analysis of the lifting case (refs. 2 to k) is complicated by 
the simultaneous appearance of lift and drag coefficients. The principal 
result of the analysis is shovn in figure 5- The vehicle is presumed to 
enter the atmosphere at some specified lift and drag coefficients indi- 
cated as point 1 on the schematic polar at the upper left. At some 
point in the trajectory, shown at the right, the resultant force is held 
constant and the coefficients are varied through the bottom of the 
pull-up (point 2) and trace a portion of the drag polar from points 1 
to 2. For this type of entry, the loading that is held during the modu- 
lation along arc 1,2 is found to vary inversely with a function T. 

(See fig. 5 -) The integral term of this function is the more important 
and is a simple line integral along the polar. Since the peak loading 
is minimum for maximum T, it is of interest to know the proper choice 
of points 1 and 2 which yields maximum T for any given polar. By 
differentiation of the r function, the maximum value of T is found 
to occur for modulation between maximum lift coefficient and minimum 
drag coefficient. 

The optimum polar with which to connect specified maximum lift and 
minimum drag points is indicated in figure 6 and is derived as follows: 
Since the endpoints are specified, it is a question of maximizing the 
integral term (fig. 6) of the function T. Since the integration is with 
respect to Cp, the largest possible at every Cp value is required. 

Lines of constant C R are circles centered cn the origin, one of which 

is indicated in the figure. The dashed curve (fig. 6) represents a 
realistic polar between points 1 and 2. For maximum at a given 

C R , the dashed curve must be pushed as far clockwise as possible along 
the lines of constant C R . This movement forces the dashed curve into 

coincidence with the right -angle polar between points 1 and 2 (drawn as 
a solid line). For the right-angle polar, T is a simple logarithmic 
function only of the ratio of maximum lift coefficient to minimum drag 
coefficient. The values of T for real polers will be somewhat less 
than the values for right -angle polars between the same maximum lift and 
minimum drag points. However, the previously mentioned conservative 
nature of the constant-velocity assumption leads in many cases to better 
loading estimates when the right-angle polar is substituted for the 
actual polar. 

For a family of vehicles with Newtonian drag polars, the results of 
the lifting analysis are shown in figure 7- The base loading is, in 
this case, that experienced in entry at the corner of the optimum polar. 
(See appendix B.) Although the analysis yields results in terms of the 


voronr 



i J H 04 ON 


7 


lift modulation parameter, in this figure the corresponding ( l / d )mAX 
of the vehicle has been used. Since vehicles of increasing (L/ d )mAX 

have increasing lift -modulation parameters, the trend of the curves is 
the same on either basis. As in the ballistic case, diminishing returns 
appear at the higher parameter values, but appreciable gains are evident 
for ( l / d )mAX values as low as unity. The ratios plotted in figure 7 

are independent of the velocity and entry angle to the order of approxi- 
mation of the analysis. 


NUMERICAL RESULTS 


Figure 8 shows some concrete examples of the acceleration levels 
to be expected for steep entry at satellite speed with and without modu- 
lation. (See appendix A.) The entry-angle range is between 6° and 12°. 

On the left are shown the unmodulated g levels for the ballistic case, 

L/D = 0, and for a pure lifting case, L/D = «>. Curves for the higher 
finite L/D values cluster in the neighborhood of the L/D = «> curve. 
With a nominal limit of lOg only the lifting vehicle at the low end of 
the range is acceptable. In the modulated case it is necessary to put 
limitations on the modulation capability. For example, in the ballistic 
case, an indefinitely large modulation capability means indefinitely 
low peak loadings. A vehicle which had and used such a capability would 
strike the earth just as if the earth had no atmosphere. Similar con- 
siderations preclude an infinite modulation ratio in the lifting case. 

For the ballistic case, then, a 50-to-l drag -modulation capability was 
assumed as a representative high value. For the lifting case an L/D 
of 2 was specified as representing a high lift capability. On the basis 
of lOg being permissible, the ballistic vehicle is acceptable through 
most of the range and the lifting vehicle is satisfactory through the 
entire range. For escape speed the absolute values of the loadings are 
different, but the relative magnitudes axe about the same. 

Until now nothing has been said about the effect of modulation on 
entry corridor width or on the heating penalties associated with modu- 
lated entry. Figure 9 deals with these considerations. (See appendix C.) 
The vertical scale is Chapman's entry corridor width in miles for a 
lOg entry at escape velocity; the horizontal scale is the (l/D) max of 

a family of vehicles with Newtonian drag polars. Two modes of operation 
are indicated by solid lines: full modulation from maximum lift coeffi- 

cient and no modulation from maximum lift-drag ratio. Rather spectacular 
gains in corridor width are evident for the higher L/D, at least in 
terms of the L/D = 0 value which is about 7 miles. Some idea of the 
relative heating penalties is indicated for a vehicle with the high L/D 
of 2.8. This represents a difficult case. Vehicles of lower (L/ d )mAX 
will have easier heating problems. 
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Trajectories corresponding to the four nodes of operation of the 
high L/D vehicle are indicated in figure 10. The dashed portions of 
the trajectories indicate the portions corresponding to the modulation 
period. Note that all the pull-ups are completed within a minute. The 
reference total heat is that absorbed down to the bottom of the pull-up 
for entry at maxim-urn lift coefficient without modulation. Figure 10 
shows this to be the highest altitude lowest heating pull-up and fig- 
ure 9 shows it to have the narrowest corridor. Entry at (l/dJj^ with- 
out modulation yields a wider corridor but 2.3 times as much heat absorp- 
t'ion. Passing to full modulation yields the widest corridor, the deepest 
penetration, and 6.1 times as much heat absorbed. The remaining trajec- 
tory corresponds to partial utilization of the vehicle capability. For 
this case more than 100 miles of corridor can be maintained for 1.7 times 
as much heat absorption as in the narrow -corridor low -heating pull-up at 
maximum lift coefficient without modulation. In terms of weight, this 
means about 70 percent more ablative material is required. In general, 
the question comes down to a trade-off between corridor width and heating 
penalties and can be settled only by a systems study of the mission. Since 
modulation from maximum lift coefficient is more favorable from both the 
heat absorption and loading standpoints, it is merely a question of how 
far to modulate. 

Results for drag modulation on a ballistic vehicle are not shown in 
figure 9 since the gains are small. At r fi = 21, calculations have shown 

the ballistic vehicle to gain in corridor width from 7 to 30 miles while 
heat absorption rises about 60 percent. The unit of heat absorption is, 
of course, different from that used in the lifting example just discussed. 

In figure 10, the peak heating rates in -the pull-ups range from 
about 500 Btu/ft 2 -sec in the pull-up at (Cl) kax without modulation 

to about 2,800 Btu/ft 2 -sec in the pull-up witi full modulation. All 
the heating values shown in figures 9 and 10 lefer to the stagnation 
point of a 1-foot -radius sphere at which point the heating rate is 
presumed to vary as p l/2 v 3.l5. 

CONCLUDING REMARKS 


In conclusion, it may be said that coefficient modulation on a 
vehicle with good lifting capability offers the possibility of sizable 
loading reductions or, alternatively, wider corridors so that steep 
entries become practical from the loading standpoint. The amount of 
steepness depends on the heating penalty that is acceptable; however, 
the price of sizable fractions of the total possible gains does not 
appear to be excessive. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Field, Va., April 11, i960. 
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APPENDIX A 

FORMULAS DERIVED BY USE OF ANALYTIC ASSUMPTIONS 


The analytic assumptions indicated in figure 3 have allowed a number 
of simple peak loading formulas to be derived. These formulas are:^ 

Unmodulated ballistic (ref. 5): 


Resultant force 
Weight 

Unmodulated lifting (ref. 2): 


Resultant force 
Weight 

Modulated ballistic (ref. l) : 


sin 7 


K 


= 2(l - cos 7 0 ) 


fW 


K 

J 


Resultant force sin 7 n „ l-A. 

— — Ae K 

Weight e 


- 1 


r B = 


(O ■g A ■g l) 


Modulated lifting (ref. 3) : 




Resultant force _ ^(l - cos 7 0 )K 
Weight p 

r = sinh _1 r L (for ideal polars) 


(Ala) 


(Alb) 


J 


The additional assumption that the net force due to gravity and coordinate 
acceleration is a small fraction of the aerodynamic force is implied by 
these formulas. Details are given in the indicated references. 


^Differences of notation from that of the references have been 
introduced. In the case of reference 2 a factor accounting for velocity 
decrement has been suppressed to conform with the constant -velocity 
assumption. 
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The quantity K is the constant acceleration in g units required 
to absorb the initial kinetic energy during straight line motion over a 

length equal to the scale height 3 \ that is. 


K = 



The formulas ( Al) above have been used to plot figure 8. In fig- 
ure 8 the parameter values assumed are: Vq = 25*900 fP s * P = 23*300 ftj 

and g = 32.2 fps. 
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APPENDIX B 

DETAILS OF CALCULATIONS OF ALLEVIATION FUNCTION 
IN LIFTING CASE 


For attitude modulation on a vehicle with an ideal polar (fig. 6) 
the minimum peak loading has been shown to be (ref. 3) 


( Resultant force 
Weight 


Mod 2 gp" 


2(1 - cos y Q ) 
r Opt 


(Bla) 


where 


Opt 


-1 


(Bib) 


The corresponding peak loading for entry without modulation (at the 
corner of the ideal polar) is (from appendix A) 


'Resultant force 
k Weight 


v 0 2 2 (! - cos y Q ) 


UniriOd 2 g 3 


-1 


(B2) 


1 + r T 


The alleviation ratio shown in figure 7 is the quotient of equations (Bl) 
and (B 2 ); that is. 


Modulated resultant force 
Unmodulated resultant force 


\jl + r L 2 sinh _1 r L 


(B 3 ) 


The Newtonian drag polars of reference 6 have been used to provide the 
relation between r^ and (L/D)^^ which is implicit in figure 7. 
These polars may be written in the form 


Cl Cd MAx( 1 " r B ) 


Cd C %ax 


'N 


sin a cos a 


— + 

(l - — |sin^a 

r B 

\ r B/ 




(B 4 a) 


J 
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The angle for maximum lift coefficient is tan~^\f2 and that for 
maximum lift-drag ratio a is defined by the relation 

(L\ = 1 + 3 cos 2a ( B 4 b) 

Vw MAX 3 sin 2a 

The corresponding ratio rg of maximum drag to minimum drag on the polar 
is found from the relation 


2 3 sin a - sin 3cc 

r B - 1 1+3 cos 2 d 

The relation between r B and r L for the po .ars (B4a) is 


(B4c) 


r T = ^S(r B - 1) (B4d) 

Eliminating r-g between equations (B4c) and (B4d) yields r^ as a 
function of a. and hence through equation (Bib), of (L/D)^^. 

Although figure 7 shows a correct trend, proper use of the numerical 
values requires that the basis of the figure oe clearly understood. The 
T values used in figure 7 are those corresponding to optimum polars 
between the maximum lift and minimum drag points of the polars given in 
eq. (B4a). Also, the unmodulated peak loading (the unit loading) is 
that experienced in an entry for which L/D = Thus, at the higher 

L/D values, the reference loading is essentially that experienced with- 
out modulation at infinite L/D. The ratios of figure 7 can be scaled 

to whatever values of r and are appropriate 

by the relations 

/ Resultant force \ a 1 
\ Weight / Mod “ r 

/ (B5) 

/ Resultant force \ ih + f— 

V Weight /Unmod V 
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APPENDIX C 


DETAILS OF CALCULATIONS OF CORRIDOR WIDTH 


The curves and points of figure 9 correspond either to numerical 
integrations or to formulas such as those used to construct figure 8. 

To avoid confusion, the construction of figure 9 requires explanation. 

The conic perigee point for the undershot boundary was determined 
in all cases by the relation: 

( r ° “ ^ undershoot = r ° Sin2(7 ° } MAX (Cl) 

which holds for parabolic orbits. The values of (7 0 )max corres P° ndin S 

to a given load limit (lOg in fig. 9) were determined either by the 
formulas of appendix A or by numerical integrations. 

The overshoot boundary was defined as the altitude at which 

in inverted flight yields 1 g unit of aerodynamic lift. This definition 
is consistent with the assumption of constant velocity. The actual over- 
shoot boundary is somewhat higher. 

The values of stagnation-point relative heat shown in figure 9 and 
the corresponding corridor widths were all obtained by numerical integra- 
tion in the ARDC atmosphere (ref. 7 ) for the L/D =2.8 vehicle of refer- 
ences 3 and A. For this vehicle (r Q - r pe ) overshoot - (0.10)l0 6 feet. 

The initial conditions were V G = 36,500 feet per second at 

h Q = 350,000 feet. The curve marked unmodulated was faired to agree 

with the machine computations for this vehicle. Results from references 6 
and 8 lie slightly below this curve. The trajectories of figure 10 were 
obtained in the same numerical integration which yielded the heating 
values. 

The curve marked "full modulation" is based on the alleviation func- 
tion indicated in figure 7 and explained in appendix B. 
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MODULATION PRINCIPLE 



Figure 1 
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Figure 2 
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ANALYTIC ASSUMPTIONS 



Figure 3 


BALLISTIC-CASE GEOMETRY MODULATION 
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Figure 4 
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LIFTING-CASE ANALYTICAL RESULTS 



Figure 5 


LIFTING-CASE OPTIMUM POLAR 
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Figure 6 
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LIFTING-CASE ATTITUDE MODULATION 


MODULATED 
RESULTANT FORCE 
UNMODULATED 
RESULTANT FORCE 



Figure 7 
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Figure 8 
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EFFECT OF MODULATION ON CORRIDOR WIDTH 



Figure 9 


PULL-UPS AT ESCAPE SPEED 

UNMODULATED MODULATED 
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Figure 10 




